On shear-free perturbations of f(R) gravity 
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Recently it was shown that il the matter congruence of a general relativistic perfect fluid flow in 
an almost FLRW universe is shear-free, then it must be either expansion or rotation-free. Here we 
generalize this result for a general f(R) theory of gravity and show there exist scenarios where this 
result can be avoided. This suggests that there are situations where linearized forth-order gravity 
shares properties with Newtonian theory not valid in General Relativity. 

PACS numbers: 04.25.Nx 



I. INTRODUCTION 

The observational evidence for the accelerated expan- 
sion of the universe, and the introduction of the concept 
of Dark Energy has put theoretical cosmology into crisis. 
This is due to the fact that despite an increasing amount 
and quality of data, no model has been proposed thus far 
that is able to give a completely satisfactory theoretical 
explanation of all the available observational data. 

Among the many different ways to achieve cosmic ac- 
celeration, the modification of the classical gravitational 
action based on General Relativity has recently gained 
much attention [3-0] • The reason for this popularity 
is due to the fact that these models provide a some- 
what more natural explanation of the cosmic accelera- 
tion: this effect is due to corrections to Einstein gravity 
which are directly related to the characteristic properties 
of the gravitational interaction. Most investigations of 
higher order gravity have focused on Fourth Order Grav- 
ity (FOG), i.e., on gravitational Lagrangians in which the 
corrections are at most of order four in the metric, and 
in what follows we will also focus on these models. 

Because the field equations resulting from FOG are 
highly nonlinear, difficult conceptual and technical issues 
arise which need to be resolved in order to uncover the de- 
tailed physics of these models. Consequently it is crucial 
to develop new methods which are able to assist in resolv- 
ing these problems. Two such approaches, the dynamical 
systems approach to cosmology and the 1+3 covariant 
approach, have proved very useful in contributing to our 
understanding of how the astrophysics and cosmology of 
these theories differ from what is found in General Rela- 
tivity. Most of the work thus far has focussed on the dy- 
namics of homogenous cosmological models Q , the linear 
growth of large scale structure |8j and on finding exact 



solutions which describe the gravitational field of stars 
and compact objects [9(. Of particular importance is to 
understand the relationship between the Newtonian and 
Relativistic limits of FOG which is important in describ- 
ing the dynamics of nonlinear fluid flows in such theo- 
ries. This is relevant both in the physics of gravitational 
collapse and the late (nonlinear) stages of structure for- 
mation (loj . 

Central to all of these problems is the differential prop- 
erties of timelike geodesies which describe the fluid flow in 
cosmology. In general the kinematics of such fluid flows 
are described by the expansion O, shear (or distortion) 
(Jab, rotation uj c 1 and acceleration A a of the four- velocity 
field u a tangent to the fluid flow lines, their governing 
equations obtained by contracting the Ricci identities 
(applied to u a ) along and orthogonal to u a , which de- 
termine how they couple to gravity [TT| . 

The delicate relationship between the kinematic quan- 
tities in Newtonian and relativistic fluid flows in General 
Relativity is most strikingly seen in a remarkable result 
first obtained by Ellis in 1967 12j. In this paper it was 
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found that if the four velocity vector field of a barotropic 
perfect fluid with vanishing pressure is shear-free, then 
either the expansion or the rotation of the fluid vanishes. 
This is a purely local result to which no correspond- 
ing Newtonian equivalent appears to hold, as counter- 
examples can be explicitly constructed It is there- 
fore interesting to ask whether such a result holds in the 
more general setting of FOG. 

As a first step towards this goal, we examine whether 
this result holds in situations where the hydrodynamic 
and gravitational equations have been linearized about 
a Friedmann-Lemaitre-Robertson- Walker (FLRW) back- 
ground Eng. These almost FLRW models can be 
thought of as lying somewhere between the full nonlinear 
situation and Newtonian theory, at least in the cosmo- 
logical context, and therefore an analysis of the theorem 
in this context could shed some light on the generality 
of the result in FOG, indeed, since it has already been 
shown in earlier work that cosmologies with a bounce oc- 
cur more naturally in such theories [2fj] , one might expect 
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a somewhat weaker version of the theorem to emerge. 

We show that if the 3-curvature vanishes, then the re- 
sult of [IH can always be avoided for vacuum universes. 
We also demonstrate there is at least one physically re- 
alistic non-vacuum case in which both rotation and ex- 
pansion are simultaneously possible. 



II. FIELD EQUATIONS FOR f(R) GRAVITY 

We know that for homogeneous and isotropic space- 
times, a "sufficiently general" fourth-order Lagrangian 
only contains powers of R and we can write the action as 



A = | / d 4 xy 



-g [f(R) + 2C T 



(1) 



where C m represents the matter contribution. Varying 
the action with respect to the metric gives the following 
field equations: 

f'G ab = T^+i(/-i?/')g ab +V b V a /'- 5ab V c V c /' , (2) 

where /' denotes the derivative of the function / w.r.t 
the Ricci scalar and T^ is the matter stress energy tensor 
defined as 



T2 = n m u a u b + p m h ab + q?u b + qTu a + tt2 



(3) 



Here u a is the direction of a timelike observer, 



u a u b is the projected metric on the 3-space 



K = i 

perpendicular to u a . Also /x m , p rn , q m and 7r^ denote 
the standard matter density, pressure, heat flux and 
anisotropic stress respectively. Equations ^ reduce to 
the standard Einstein field equations when f(R) = R. 

We can write the modified field equations ([2]) as 

Gab = T a l +T a R b = T ab , (4) 

where 

rpm 

/' 



rpr) 



(5) 



Tl = jj[Hf- R f)9ab + V b V a /' - ,g a6 V c V c /'] , (6) 

and the thermodynamic quantities of this fictitious 11 cur- 
vature fluid' are given by 



R _ T R a b R _ luabrpR 
P — 1 ab u U i P — 3 n 1 ab > 



n R _ ubrpR c R _ 
la — a a 1 bc U I 71 ab ~ 



d — \h ab h c d 



rpR 

1 cd 



(7) 



(8) 



The total thermodynamic quantities are then written as 



fx'" 

T 



p" 1 



Q 



T^ab 



r 



'ab- 



(9) 



Using the standard 1+3 covariant approach [LT| . two 
derivatives are defined: the vector u a is used to define 
the covariant time derivative (denoted by a dot) for any 
tensor T°" c ..d along the observers' worldlines: 



J- c.d — U V e J c ..d 



(10) 



and the tensor h a b is used to define the fully orthogonally 
projected covariant derivative V for any tensor T a ^ ' c ..d'- 

VeT a - b c .. d = h a f hP c ...h b g hi d h r e V r Tf-\.. q , (11) 

with total projection on all the free indices. Angle brack- 
ets denote orthogonal projections of vectors and the or- 
thogonally projected symmetric trace-free PSTF part of 
tensors: 



y<«> = h a b V b , T< ah > = hl a c h b) d - \h ab h cd 



T* 



cd 



(12) 



This splitting of spacetime also naturally defines the 3- 
volume element 

eabc = -V\9\Sf a S 1 b 6 2 c 6%u d ) (13) 

with the following identities: 

e abc e def = 3% h e b h f c] , e abc e dec = 2\hf a h e b] . (14) 

The covariant derivative of the timelike vector u a can 
now be decomposed into the irreducible parts as 



V a it 6 = -A a u b + hh ab Q + a ab + e a . 



(15) 



where A a — it a is the acceleration, 9 = V a it a is the 
expansion, a a b — V uU b \ is the shear tensor and u> a = 
e abc \7 b u c is the vorticity vector. Similarly the Weyl cur- 
vature tensor can be decomposed irreducibly into the 
Gravito-Electric and Gravito-Magnetic parts as 



E a .b = C abcd u c u a 



~ J (ab) , H a b — \tacdC be U e — £f( ab ) , 

(16) 

giving a covariant description of Tidal forces and Gravi- 
tational radiation respectively. 



III. LINEARIZED FIELD EQUATIONS ABOUT 
FLRW BACKGROUND 



In the 1+3 covariant perturbation theory [14H19L 1 211 . 
HH, the quantities that vanish in the background space- 
time are considered to be first order and are automati- 
cally gauge- invariant by virtue of the Stewart and Walker 
lemma }23| . We consider the background to be FLRW 
where the Hubble scale sets the characteristic scale of 
the perturbations. In the perturbed spacetime the stan- 
dard matter is considered to be a perfect fluid with the 
Energy Momentum tensor given by: 



T% = (» m +p m Wu b +p m 9ab. 



(17) 



3 



Furthermore, we assume standard matter to have a 
barotropic linear equation of state p m = wfi m satisfy- 
ing the Weak and Dominant energy conditions: 



H m > ; fi m +p m > ; M m > \p m \ 



(18) 



Since the matter is a perfect fluid the heat flux (g™) 
and the anisotropic stress (7r^) vanish in the perturbed 
spacetimc. In addition, since we consider shear-free per- 
turbations the shear tensor a a t vanishes identically. 

For the Curvature Fluid the linearized thermodynamic 
quantities are given as 



// 



1 

J' 



URf'-f)-ef"R + f"V 2 R 



(19) 



Constraint equations 



(C ) ab := E ab - V {a A^ - \nf = , (29) 



(d) a := V a 6 - h abc V b uj c - \q a R = , (30) 



(C 2 ) := V a uj a = 



(C 3 ) ab := H ab + V {a uj b '> = . 



(31) 



(32) 



p* = ± 



AS -R.n + f"R + f'"R 2 



(ef'R - ,/"v 2 i?) 



(20) 



(CiT :=V> m + ( Mm +p m )A a = 0, (33) 



(C 5 ) a := V b E ab + iV fc 7T^ - ±V> + ±6<& = , (34) 



q B = -- 



f"RS7 a R + f"V a R-\f"@V a R , (21) 



(22) 



With the conditions above, the linearized field equations 
are then given by: 



Propagation equations 



e-v a ^ a = -ie 2 -i(/i + 3p) , (23) 



Ja) _ i e abc V b A c = -|6w a , 



E (ab) _ ^(ay^b = _Q E ab _ l^ab 
~2 V 1R ~ 6 B7r fl ' 



(24) 



(25) 



H (ab) + € cd{ay cE )b = _Q H ab + l £ cd<a V c ^ 6 , (26) 



(C 6 ) a := V 6 ff ah + (/x + P )u a + \e abc V b q* = . (35) 

We note that the constraints (Ci) a , (C 2 ), (C 3 ) afc , (C 5 ) a 
and (Ce) a are the constraints of Einstein field equations 
for general matter motion and are shown to be consis- 
tently time propagated along u a locally in General Rela- 
tivity. However, the conditions a ab — and = give 
the two new constraints (Co) ab and (64)° respectively. 
In what follows we will use the following linearized com- 
mutation relations for shear-free congruences. For any 
scalar 4> 

[V a V 6 - V 6 V a ]</> = 2e afcc a; c 0, 

e abc V b V c <t> = 2oj a <p. (36) 
If the gradient of the scalar is of the first order, we have 
[V a V 6 V a - V 6 V 2 ]0 = \RV b 4> (37) 

and 

[V 2 V 6 - V 6 V 2 ]0 = ii?V b + 2e dbc V d (^), (38) 

where R = 2 (/i — |0 2 ) is the 3-curvature scalar. Also 
for any first order 3- vector V a — V^ a \ we have 



[V a V b -V b V a ]V a = \Rh a [a V b 



(39) 



(Mm +Pm)6, 



(27) 



(40) 



A + v a c = -(/i + p)e ; 



(28) 



h a c (v 2 v c )- = v 6 (v ( V a> )- - ±ev 2 y a . 



(41) 
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IV. CONSISTENCY OF THE NEW 
CONSTRAINTS 



Again using ([36]) and (fT4]) in the above equation and 
linearizing we get 



We have already seen that the conditions of shear-free 
perturbations together with the matter being described 
by a perfect fluid in the perturbed spacetime, gives the 
new constraints (Co) ab and (64)° respectively. To check 
their compatibility with the existing constraints of Ein- 
stein's field equations, we substitute (Co) bd into (Cs)b to 
obtain 



W d W {b A d) - i V b /i + V d 7r* + §0g* = 
Now from the constraint {C^) b we have 



10 + 1 



(42) 



(43) 



where = ln(/i m ). Using equation (j43|) in (|42|) we get 
the constraint 

V d V (b V d> + § V bM - VVd - S0?6 R = 0. (44) 



w + 1 



We note that for the new constraints to be compatible 
with the existing ones the above constraint must be sat- 
isfied. To check the spatial consistency of the above con- 
straint on any initial hypersurface we take the curl of (|44[) 
to get 



-e V C V ir bd - jHe V c q b = , 
which, on using (j3"o]l gives 



(45) 



10 + 1 



-e^VrVVrtV" 



-±e e QCb v c 



e ac6y yd 



y-RVbR+jfVbR 



f"~ ■ e/"~ 



3/' 



V b i? 



(46) 



Breaking the PSTF part according to equation (jl2|) . us- 
ing the commutators (|37]l. (|38]l and keeping only terms 
up to first order, we have: 



10+1 



§ V c V b V 2 + ±i?V c V b + (^fcVcV'V 



+| W Q /i - ^-e acb V c V d V (b V d) i? + i-_ e - b V c V b ii 



3 -e acb V c V b i? - L— e acb V c V b R = 0. 



3/' 



9/' 



(47) 



w + 1 



fifcj°<£ - 0V fc V Q w fc + 0V 2 w° 



+ |w a /i 



f" 



§ V c V b V 2 i? + ±RV c V b R + Re dbk V c V d w k 



f 

This can also be written as 



f'OR 2 - \f"Q 2 R + f'GR 



(48) 



w 



w + 1 



2 .Rw a </> - (jN k V a u k + 0V 2 w 



72, ,a 



-i? 



/' 

2w a 
3J 7 



+ |w a /i 



/'"9i? 2 - i/"6 2 i? + /"9i? = 0. (49) 



Now, from relation (|3"5|) and using (j3"Tj) we know that 

V k V a co k = \Ruj a , (50) 



and from (|27p we have 

= -(1 + 10)6 
Thus rearranging terms gives 
R 



(51) 



we 



^ +P )eiu a + R^ r 



+V 2 cj a 



29cj a 



f'R 2 - i/"6i? + f'R 



I- 

= . 



(52) 



Since 



+ p 



/' 



3/' /' 



/'"i? 2 /" 



/' 3/' 

the above equation, to linear order, simplifies to 

2(1 + w)n m Quj a 



(53) 



we 



~W° + v 2 w a 
3 



+ 



3/' 



+R 



r 



\Rw a + V 2 w a 



(54) 



Further manipulation leads to 



r^T \ 5 2(i + z0)^ m e 

3 3/' j + 3/' 



Rf" 
f 



wQ V 2 u a = 0. 



(55) 
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We know that in terms of the scale factor a(t) of a FLRW 
spacetime, the expansion, acceleration, jerk and snap pa- 
rameters are defined by the following relations: 



where we have used 



e = 3- 



aa 



aa 

q = -to-) 

„3 J4 



J 



d 4 a 
U~dF' 



(56) 



(57) 



From the above equations we can easily see that the time 
propagations of these quantities can be written as 



e 



i 



i 



e 2 (i + 9 ) , 



q = -~®{j-q-2q 2 ) 

e = ie 3 (2 + 3<z + j) , 



j = -e(s + 2j + 3qj) , 

q = -\® 2 [s + 2j - 3q 2 + 6qj - 6q 3 ] 

Then, the Ricci scalar R is given by 
2 



and hence 



R=^Q 2 (l-q)+R 



where 



Q = \@ 2 {j-q-2) + k, 



(58) 
(59) 
(60) 
(61) 
(62) 

(63) 
(64) 
(65) 



(4 + 5q + j + jq + s)9 2 + 6R . (66) 
This means that we can rewrite (1551) as 



r 2 



V 2 uj a 







(67) 



We can see from this equation that spatial consistency 
requires the vanishing of either O or the terms in the 
curly brackets. 

To check for temporal consistency of the new con- 
straint (|67l) we take its time evolution which can be writ- 
ten as 



io a { 6 



(1 - w)P - (l + w) (3w + 5)/' + 4/"Q 

■■K H - -nr. fJ-r, 



r 



6/' 



P 
P 



l + w 

(— )Mr, 



= 0, 



(68) 



w a = (w - |)ew a 



and 



(V 2 ^)- = (3w „ 5) 9VV 
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We have also defined 

P : 



6 



eRu a 



f" 

From (|66[) . we can write 



3w 
~2~ 



(69) 



(70) 



(71) 



(72) 



where 



2 // 



z = f(T- ( 7 )2 J 03 

f" / 

+ 9f 7 (( 4 + 5 9 + J+J9 + s )© 2 



6i? 



(73) 



Equation (|67[) can then be rewritten as 

(1 + w) (3w + 5)/' + 4/'' 



(1 - w)P ~ 



f 



6/' 



J Mm 







(74) 



It follows that for the new constraints to be spatially 
and temporally consistent we must have either w a O = 
or the expression in the curly brackets must vanish. It 
is interesting to see whether there exist solutions of a 
given f(R) theory of gravity which can avoid the Ellis 
condition. 

From (fT4")) . it is easy to see that if the 3-curvature van- 
ishes, then the result of [24[ can always be avoided for 
vacuum universes (/i TO = 0). This implies that a shear- 
free, spatially flat vacuum universe in any f(R) theory 
can rotate and expand simultaneously in the linearized 
regime. 

The non-vacuum case is more difficult to analyze in 
general; however, as we will see below there does exist 
at least one non-trivial case which does violate the Ellis 
condition. 

For a flat Milne universe, where the matter energy den- 
sity is given by (j. m = 3 ffi m) , we have 



9 
R 



±e 2 

3 ' 



|e 2 



a(R) = 



R 



R 



2 i?2 . 



(75) 
(76) 

(77) 
(78) 
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Substituting these quantities into the Friedmann equa- 
tion 



ie 2 = l 

3 /' 



f-l>m 



Rf'-f 



QRf" 



one gets 



R 



2 d 2 f(R) , f(R) 



Mo 



dR 2 2 a(i?) 3 ( 1 + u ') " 

which has the following general solution: 

4^0 



(79) 



(80) 



f(R) = C.R 1 ^ + C 2 R^ 



1 + 12w + 9w 2 



R— 



(81) 

Let us only consider the particular solution (the last term 
of the above equation) , which is an R n - theory of gravity. 
Now, if we look at (l74l) . for the corresponding flat Milne 
universe in R n gravity, the term in the curly brackets 
reduces to 



(1 + w)Hr, 

6/' 



[3w + 9 - An] = 0. 



(82) 



Comparing solutions (|8"2"j) and the particular solution of 
dHJ) (with n = 3(1 + w)/2) we get w = 1 if \i m ^ 0. In 



other words, for a stiff fluid in R 3 gravity, there exists 
a flat Milne-universe solution which can rotate and ex- 
pand simultaneously at the level of linearized perturbation 
theory. 



DISCUSSION AND CONCLUSION 



In this paper we consider shear-free fluid flows in f(R) 
gravity in situations where the hydrodynamic and grav- 
itational equations have been linearized about a FLRW 
background. This extends recent work by Nzioki et. al, 
which considered such situations in General Relativity 
24j . We showed that if the 3-curvature vanishes, then 
the result of [IH, [25| can always be avoided for vacuum 
universes. We also demonstrated there is at least one 
physically realistic non-vacuum case in which both ro- 
tation and expansion is simultaneously possible. This 
suggests that there are situations where linearized forth- 
order gravity shares properties with Newtonian theory 
not valid in General Relativity. 
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